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Multigap superconductivity and strong electron-boson coupling in Fe-based 
superconductors: A point-contact Andreev-reflection study of Ba(Fex x Co x ) 2 As 2 single 

crystals 
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J. Jiang, 3 J.D. Weiss, 3 E.E. Hellstrom, 3 and R.S. Gonnelli 1 
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Directional point-contact Andreev-reflection (PCAR) measurements in Ba(Fei_ x Co x )2As2 single 
crystals (T c =24.5 K) indicate the presence of two superconducting gaps with no line nodes on the 
Fermi surface. The PCAR spectra also feature additional structures related to the electron-boson 
interaction, from which the characteristic boson energy Qb(T) is obtained, very similar to the spin- 
resonance energy observed in neutron scattering experiments. Both the gaps and the additional 
structures can be reproduced within a three-band s± Eliashberg model by using an electron-boson 
spectral function peaked at Qo = 12 meV ~ fib(O). 

PACS numbers: 74.50.+r , 74.70.Dd, 74.45. +c 



The discovery of the first class of non-cuprate, Fe- 
based high-temperature superconductors in 2008 brought 
great excitement in the scientific community [lj. The 
phase diagram of these compounds (although still im- 
perfectly known) looks similar to that of copper-oxide 
superconductors Q and, as in cuprates, superconductiv- 
ity emerges "in the vicinity" of a magnetic parent com- 
pound. The electron-phonon interaction seems not to be 
sufficient 3 to explain their high T c (up to 55 K [J) 
even by considering a magnetic ground state Q . A spin- 
fluctuation-mediated pairing mechanism has been early 
proposed instead, which predicts the occurrence of a sign 
change of the order parameter on different sheets of the 
Fermi surface (s±-symmetry) 0]. This picture is nat- 
urally based on the proximity of the superconducting 
phase to a magnetic one, on the existence of disconnected 
Fermi surface (FS) sheets, and on the multiband charac- 
ter of superconductivity in these compounds, which are 
nowadays almost universally accepted [2]. The s± model 
itself is strongly supported by various experimental re- 
sults § which indicate the existence of multiple nodeless 
gaps on different sheets of the FS, although the possible 
emergence of gap nodes in some systems, along certain 
directions or in particular conditions d, M is still de- 
bated. The role of spin fluctuations (SF) in the pairing 
has also found support in neutron scattering experiments 
that have revealed a spin resonance energy which scales 
linearly with T c ,2\. Finally, it has been recently shown 
that a multiband s± Eliashberg model can reproduce sev- 
eral experimental quantities (such as gaps, T e , kinks in 
the band dispersion and effective masses [ll|, [l2() by as- 
suming that the mediating boson has a characteristic en- 
ergy similar to the spin-resonance one. 
In this paper we report on directional PCAR measure- 
ments on high-quality single crystals of the e-doped 122 



compound BaFe1.sCoo.2As2. The results prove the ex- 
istence of two superconducting gaps with no line nodes 
on the FS, and whose amplitude is almost the same in 
the ab plane or along the c axis. The PCAR spectra 
also present structures that can be related to a strong 
electron-boson interaction (EBI). The characteristic en- 
ergy Hi, of the mediating boson extracted from the PCAR 
curves decreases with temperature and is very similar to 
the resonance energy of the spin excitation spectrum [i~3| . 
Moreover, both the gaps and the additional EBI struc- 
tures in the PCAR spectra can be reproduced within an 
effective three-band s± wave Eliashberg model using a 
boson energy Slo = 12 meV ~ O&(0). All these results 
strongly support a spin-fluctuation-mediated mechanism 
for superconductivity in this compound. 
The BaFe1.sCoo.2As2 (10% Co) single crystals were pre- 
pared by the self-flux method [14| under a pressure of 280 
MPa at the National High Magnetic Field Laboratory in 
Tallahassee. The typical crystal sizes are w 1 x 1 x 0.1 
mm 3 . The onset of the resistive transition is T° n = 24.5 
K with AT C (10%-90%) = 1 K (see inset to FigHJ. In- 
stead of using the standard technique where a sharp 
metallic tip is pressed against the material under study, 
the point contacts were made by putting a small drop of 
Ag paste on a fresh surface exposed by breaking the crys- 
tal. Contacts made in this way are very stable and the 
differential conductance curves, obtained by numerical 
differentiation of the I-V characteristics, can be recorded 
up to « 200 K [l5| . As an example, Fig Q] shows the 
raw conductance curves, recorded up to 180 K, of a 
Ag/BaFei.8Coo.2As2 point contact (i?N = 25 O) with 
current injection along the c axis ("c-axis contact"). The 
clear signatures of AR in the low-T curve and the absence 
of heating effects or dips [lj| indicate ballistic conduc- 
tion through the point contact, so that energy-resolved 
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FIG. 1: (color online) Temperature dependence of the dif- 
ferential conductance curves in a Ag/BaFei.8Coo.2As2 c-axis 
point contact. The curves are vertically offset for clarity. The 
insets show the superconducting transition as seen by AC 
magnetic susceptibility (left) and DC resistance (right) mea- 
surements. 

spectroscopy is possible. A closer inspection reveals that 
the maxima in the low-T curves present fine structures 
(indicated by arrows in FigUJ suggesting multiple gaps. 
The Andreev signal decreases on increasing T and com- 
pletely disappears at the critical temperature of the con- 
tact, T C A = 22.6 ± 0.2 K, leaving a slightly V-shaped 
normal state. On further heating, the normal-state curve 
progressively fills and completely flattens at « 140 K, the 
temperature where the long-range magnetic order sets in 
in the parent compound. Similar behavior was observed 



in 1111 Fe-based superconductors [15l.ll7|. 

In order to compare the experimental curves to a suit- 
able model, all the raw conductance curves at T < T C A 
were normalized by the normal-state curve at T C A . Fur- 
thermore, to get rid of the well-known asymmetry of 
the PCAR spectra of Fe-based compounds [HI, G3, EH 
the normalized conductance G(V) was symmetrized, i.e. 
G{V) = [G cxp (V) + G exp (-F)]/2. This preserves and 
enhances the structures we are interested in (gaps and 
EBI). The asymmetry of the original curve was taken 
into account as a source of uncertainty on the gap val- 
ues. The resulting conductance curves were then fit- 
ted to a two-band BTK model [29[ taking into account 
broadening effects and the angular distribution of the 
injected current [lfij . In this model the normalized con- 
ductance is the weighed sum of two BTK terms G(V) = 
w\G\{V) + (1 — wi)G2(V), where wi is the weight of con- 
tribution 1. Each term G\ is described by a gap value A;, 
a broadening parameter Tj (here mostly due to inelastic 
scattering in the vicinity of the contact) and the parame- 
ter Z\ which accounts for the height of the barrier at the 
N/S interface and the Fermi velocity mismatch [lill ]. 

Fig. [2] shows the setup for PCAR measurements with 
current injection along the c axis (a) and along the ab 
plane (d). Examples of normalized conductance curves 
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FIG. 2: (color online) (a,d) Sketch of c-axis and afe-plane 
contacts. (b,c): normalized conductance curves at 4.2 K for 
c-axis contacts (symbols) and their two-band fit (solid lines) 
with the relevant gap values Ai and A2. Arrows mark the 
structures related to the gaps and to the EBI. (e,f): the same 
for two afe-plane contacts. In (b) and (e), a single-band fit is 
also shown (dashed lines) with the relevant gap amplitude A. 

at 4.2 K are shown in (b) and (c) for c-axis contacts 
and in (e) and (f) for a&-plane contacts. All the PCAR 
spectra show peaks at w 4 meV and shoulders at « 9- 
10 meV. Additional structures are reproducibly present 
at 18-20 mV, although more pronounced when the An- 
dreev signal is higher. In few cases (panel f) they are 
masked by small dips, which however do not affect the 
very clear two-gap structures at lower energy. Fig. [5] 
(b) and (e) (as well as the inset to Fig(3]) clearly show 
that a one-gap BTK model (dashed line) is unsuited to 
reproduce the experimental data while a two-gap model 
allows a good fit of the experimental curves, apart from 
the structures around 20 mV. The resulting amplitudes 
of the gaps Ai and A2 are indicated in the labels. In 
all the two-gap fits of this paper w\ = 0.5 ± 0.1 and, at 
low T, r/A = 0.5 — 0.7. Finally, Z and w\ are constant 
with temperature while T is almost constant or slightly 
increases with T 15, lj|. 

From the two-gap fits of various curves we obtained the 
average values: A^ = 4.1 ± 0.4 meV and A^ = 9.2 ± 1.0 
meV for c-axis spectra and A^ b = 4.4 ± 0.6 meV and 
Ajf 3 = 9.9 ± 1.2 meV for a6-plane contacts. These re- 
sults can be compared to ARPES experiments [l9[ , which 
show two nodeless gaps in the k x k y plane. The small 
gap, located on one of the electron FS sheets, is in very 
good agreement with our Ai. Our value of A2 is instead 
about 30% bigger than the large ARPES gap, located 
on the hole FS sheet. The reason of this discrepancy 
will become clear in the following. In this concern, note 
that, although directional PCAR measurements are not 
fc-resolved, they allow probing the gaps also along the k z 
direction, not easily accessible to ARPES measurements. 
The absence of zero-bias conductance peaks (ZBCP) 
along either direction in the PCAR spectra rules out line 
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nodes on the FS both along the c axis and in the ab 
planes, but does not exclude deep gap minima or even 
zeros in small regions of the Brillouin zone 2H 21 1- 
The fact that w\ is almost independent of the direction 
suggests an almost equal degree of three-dimensionality 
of the various FS sheets in Ba(Fei_ x Co x )2As2, as also 
shown by ARPES [Hj], X-ray Compton scattering 23 1 
and first-principle calculations 0, [24| . 
Fig. [3] shows the temperature dependence of the normal- 
ized conductance of Fig. [2jc) (symbols) and the relevant 
two-band BTK fit (lines). The two-band model fits very 
well the PCAR spectrum at low T (see left inset) giving 
Ai(0) = 3.8 meV and A2(0) = 8.2 meV, which correspond 
to 2Ai/fc B T c «3.9 and 2A 2 /fc B T c w 8.5, both above the 
BCS weak coupling ratio. The temperature dependence 
of the gaps is shown in the right inset (symbols). 
It has been recently shown that in La-1111, Sm-ffff 
and Bai_ x K x Fe2As2 the experimental gap values and 
their temperature dependence can be reproduced within 
a three-band s± Eliashberg model [ll], 12 1, while two- 
or three-band weak-coupling BCS models cannot do the 
same. In Ba(Feo.9Coo.i)2As2 we can simplify the elec- 
tronic structure, according to ARPES measurements [lj| , 
by taking one effective hole band (band 1) and two elec- 
tron ones (band 2 and 3, corresponding to the outer and 
inner electron barrels in the FS as defined in Ref.[20]). 
We disregard the small hole pocket at T, predicted by cal- 
culations but not observed by ARPES. Phonons mainly 
provide intraband coupling but their contribution is ex- 
pected to be small [J while spin fluctuations (SF) 
mainly provide the interband coupling. We thus set 
AP ; h = 0.2 [31 and Ajf = A?f = so that the electron- 
boson coupling matrix becomes: 

A ph A12 A13 
Ai2^i2 A? h 
Ai 3 ^i3 A? h 

where u 12 = N 1 (0)/N 2 (0), v lz = JVi(0)/iVa(0). N t {Q) is 
the normal density of states (DOS) at the Fermi level 
for the ith-band, calculated from the first-principle LDA 
bands of the 8% Co-doped compound (25|, first shifted 
downward in energy and then renormalized by a factor 
2 to agree with the ARPES results [H, HjJ. To satisfy 
the conservation of the total charge, the energy shift is 
30 meV for the h-bands and 46 meV for the e-bands. 
Finally, the total DOS of electron bands is divided in a 
4:1 proportion between bands 2 and 3. This is consistent 
with the Raman data 2l| that suggest the existence of 
"hot spots" (where the gap is substantially suppressed) 
which occupy, crudely spea king , about 1/2 or less of one 
out of two electron pockets [20( . This uneven splitting of 
the DOS is very important to obtain a satisfactory agree- 
ment between the experimental data and the results of 
the Eliashberg model. Following the above, v\2 = 1-12 
and V13 = 4.50. As for the electron-SF spectral function, 
we used a Lorentzian curve peaked at Oij = 17q = 12 meV, 
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FIG. 3: (color online) Temperature dependence of the normal- 
ized conductance of Fig^c) (symbols) and the relevant two- 
band BTK fits (lines). All curves except the bottom one are 
vertically offset for clarity. The corresponding gaps are shown 
in the right inset (symbols) compared with the BCS-like tem- 
perature dependencies. Left inset: zoom of the curve at 4.2 K 
(symbols) with two possible one-gap BTK fits (dashed lines) 
and the best two-band BTK fit (solid line). 

in agreement with neutron scattering experiments Q. 

The only two free parameters of the model are A12 and 
A13 which are chosen so as to reproduce the experimen- 
tal gaps as well as possible [11| . The obtained gap values 
are Ai = 6.1 meV, A 2 = -3.8 meV and A 3 = -8.0 
meV (with a theoretical T c w 29.7 K). Ai (hole FS) 
and A 2 (outer electron FS) are in very good agreement 
with the ARPES experiments which actually mea- 
sured the gap only on one of the two electron FS sheets. 
Also, A 2 and A3 are consistent with the gap values ob- 
served in our PCAR experiments; resolving the inter- 
mediate gap by PCAR is a challenging task. Thus, the 
whole set of data from ARPES, PCAR and calculations 
looks consistent. The coupling constants are A12 = 0.61 
and A13 = 1.22 corresponding to a total effective cou- 
pling constant A c ff = 1.93, which indicates, as expected, 
a strong-coupling character for this compound. 

Let us now discuss the aforementioned additional struc- 
tures at about 20 mV that are reproducibly observed 
in the PCAR spectra (see figH]), and that disappear at 
the critical temperature of the contacts. We will show 
here that these structures are the signature of the strong 
electron-boson coupling, where the boson characteristic 
energy is the spin-resonance energy observed by neutron 
scattering. Figure [3](a) shows the normalized conduc- 
tance at 4.2 K of a afe-plane contact where the AR signal 
is particularly high (« 30%), and the structures at ~ 20 
mV are clearer than usual, which makes this curve partic- 
ularly interesting for our discussion. The solid line is the 
theoretical PCAR spectrum obtained from a three-band 
BTK model by replacing the constant BCS gaps with the 
energy-dependent gap functions (for details on this pro- 
cedure see § 4.3.5 of Ref. M) calculated within the same 
Eliashberg model and with the same parameters dis- 
cussed above. In the absence of a theoretical way to ac- 
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FIG. 4: (color online) (a) Comparison between an experi- 
mental AR spectrum (symbols) and the theoretical one (line) 
obtained from Eliashberg and BTK calculations (see text), 
(b) Experimental (full symbols) and theoretical (solid line) 
-d 2 I/dV 2 vs. V curves obtained from the data in panel (a). 
Open symbols: the —d 2 I/dV 2 curve averaged over 5 contacts. 
Dashed line: the electron-boson spectral function (shifted in 
energy by » A max ) used in the three-band Eliashberg calcu- 
lations, (c) Evolution of the —d 2 I/dV 2 vs. V curves with 
temperature showing the displacement of the bosonic struc- 
tures. The energy of the peak E P (T) and the corresponding 
characteristic boson energy Qb(T) axe shown in panel (d). 
Lines are only guides to the eye. 

count for the broadening parameter T within the Eliash- 
berg theory, the diffusive normal metal/superconductor 
junction model was used to adjust the amplitude of the 
curve to the experimental one 27| without changing the 
position or shape of its features. This requires fixing a 
single parameter Rd/Rb = 1-015 where Rd (Rb) is the 
resistance of the diffusive bank (of the junction). The 
theoretical AR spectrum clearly shows high-energy struc- 
tures very similar, in position and in amplitude, to the 
experimental ones. 

Fig. IDJb) reports the —d 2 I/dV 2 curve for the experimen- 
tal (full symbols) and theoretical (solid line) conductance 
curves shown in figure HJa). In low-transparency (large 
Z) point contacts on strong-coupling superconductors, 
peaks in —d 2 I/dV 2 correspond to peaks in the electron- 
boson spectral function. In the case of small Z, a small 
relative shift is observed [l6[, but here it turns out to 
be negligible (< 0.2 meV). A peak in the experimental 
—d 2 I/dV 2 is clearly visible at about 21 meV (and is ob- 
served also in the theoretical curve). Other structures 
appear around 27 mV and 40 mV. All these structures 
exist also in the —d 2 I/dV 2 curve obtained by averag- 
ing over 5 different contacts (open symbols). The energy 
of the first maximum, E p , agrees well with the energy 



of the peak in the Lorentzian electron-boson spectrum 
used in our calculations, shifted by ~ A max (dashed line) 
further indicating that a bosonic mode at £Iq is re- 
ally playing a major role in the coupling. The structures 
at higher voltage that do not appear in the theoretical 
—d 2 I/dV 2 (solid line in Fig. Hp) may be due to the ac- 
tual shape of the electron-SF spectral function and/or to 
non-linear strong-coupling effects. FigHJc) shows that, 
on increasing temperature, all the EBI structures shift 
to lower energy. Figj4jd) reports the maximum and min- 
imum values of E p over the different —d 2 I/dV 2 curves 
(full symbols) , and of the quantity E p — A max (open sym- 
bols) as a function of temperature. Note that the latter is 
the energy of the "resonant mode" in the electron-boson 
spectrum, Qb (fib ~ fig at low T) ^28|j and its behavior is 
indeed very similar to that of the spin resonance energy 
measured by neutron scattering experiments 13 1. 
In conclusion, we have shown that PCAR measurements 
give direct and clear evidence for multiband strong cou- 
pling superconductivity in Ba(Fei_ x Co x )2As2. They also 
allow extracting the characteristic energy of the mediat- 
ing boson and its T dependence, that both coincide with 
those of the spin resonance measured by neutron scatter- 
ing experiments [pjj]. This brings unambiguous evidence 
for a spin-fluctuation-mediated s± mechanism of super- 
conductivity in this compound. 
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sions. This work was done within the PRIN project No. 
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